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Abstract 

We generate a complete basis of magnetic moment operators for the N c = 3 ground-state baryons 
in the 1/N C expansion, and compute and tabulate all associated matrix elements. We then compare 
to previous results derived in the literature and predict additional relations among baryon magnetic 
moments holding to subleading order in 1/N C and flavor SU(3) breaking. Finally, we predict all 
unknown diagonal and transition magnetic moments to < 0.15 /i n accuracy, and suggest possible 
experimental measurements to improve the analysis even further. 
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I. INTRODUCTION 



The generalization of quantum chromodynamics from 3 to N c > 3 color charges, called 
large N c QCD, has opened a path to substantial progress in understanding strong interac- 
tions at both the formal and phenomenological levels. Formal successes spring from the fact 
that large N c QCD exhibits a well-defined limit, meaning that the renormalization group 
equations remain finite and nontrivial as N c — > oo. Furthermore, the counting of explicit 
N c factors organizes QCD Feynman diagrams into topological classes of decreasing signifi- 
cance with increasing powers of 1/N C , which defines the 1/N C expansion. Phenomenological 
successes build on these formal 1/N C power-counting results, but add one crucial extra ingre- 
dient: Observables calculated to appear at 0(1/N C ) or 0(1/ N^) are empirically found to be 
a factor 3 or 9, respectively, smaller than corresponding quantities calculated to be 0(N®); 
this means that even QCD with N c as small as 3 belongs to the class of theories for which 
the 1/N C expansion is meaningful. We note only that the literature to date that provides 
evidence substantiating these statements has become so extensive, that nothing short of a 
review article [l| can do it justice. 

Nevertheless, a multitude of problems utilizing the 1/N C expansion, even for well-known 
observables, remain unsolved. In this paper we focus on one very specific such set, the 
magnetic moments of the u, d, s baryons in the ground-state multiplet. In the case of large 
N c , this multiplet consists of a tower of states [2] completely symmetric under combined spin 
and flavor transformations, thus providing justification for the group-theoretical aspects 
of the old three-flavor SU(6) classification for baryons. The nonstrange members of the 
multiplets in this tower carry the (/, J) quantum numbers (|, |), (|, |), . . ., (4^, ^). The 
first (/, J) multiplet represents nucleons, which reside in an SU(3) multiplet that is an octet 
for iV c = 3; the second represents A resonances in an SU(3) multiplet that is a decuplet 
for N c = 3. Here we continue to use the SU(3) labels 8 and 10, despite the fact that the 
corresponding SU(3) representations are much larger for N c > 3 [3]. Of course, for N c = 3 
this tower truncates after the A's. While the mass of each baryon is O(N^), mass splittings 
between two low-lying states in the tower [i.e., I = J = 0(N®)] is 0(1/ N c ) |4j, supporting 
the notion of a true degenerate spin-flavor multiplet. In fact, it is only because our universe 
is somewhat closer to the chiral limit than the large N c limit that A and its partners in the 
SU(3) decuplet are unstable under strong decays: m n = 0(m u> d) = 0(N®) < m^ — rriN = 



1 



0(1/N C ). 

In a complete analysis organized according to 1/N C , the whole set of states in the spin-| 
8 and spin-| 10 [and SU(3) multiplets associated with spin |, |, . . ., =rf , which decouple for 
N c = 3], must be considered together as a single, completely symmetric spin-flavor multiplet 
with N c fundamental representation (quark) indices; we continue to denote this multiplet by 
the old SU(6) label 56, although again for N c > 3 the dimension of this representation is much 
greater. The instability of spin-3/2 baryons is taken into account simply by maintaining finite 
values for and 1/N C in the full Hamiltonian. 

We hasten to add that magnetic moments for baryons in the 56 have been considered 
in the 1/N C expansion in the past — in fact, in papers dating back a decade or more. There 
are three papers in particular that have examined these magnetic moments in the 1/N C 
expansion: Jenkins and Manohar (JM) j^, Luty, March- Russell, and White (LMRW) 
and Dai, Dashen, Jenkins, and Manohar (DDJM) |7j. Each of these papers contains a 
scheme for including a particular set of operators that contribute to magnetic moments, and 
each is discussed in detail below, once we establish a notation to describe the formalism. 

In short, however, the essential improvement of the current work over previous papers 
is completeness. Once all relevant baryon states are assumed to fill a complete spin-flavor 
multiplet — in this case, the 56 — then only a finite number of operators exist with distinct 
spin-flavor transformation properties that can generate nonvanishing baryon bilinears in a 
Hamiltonian. This number precisely equals the number of distinct observables associated 
with the given quantum numbers. For example, the 56 allows precisely 19 linearly in- 
dependent mass operators (those with AJ = 0, AJ 3 = 0, AY = 0, AI 3 = 0, T even) with 
distinct spin-flavor properties, corresponding to the masses of the eight 8 baryons, the ten 10 
baryons, and the spin-singlet S°A mixing. In the magnetic moment case (AJ=1, AJ 3 = 0, 
Ay = 0, AJ 3 = 0, T odd), one finds 27 linearly independent operators, corresponding to 
the eight 8 baryons, ten 10 baryons, the AJ = 1 S°A mixing, and eight SU(3)-breaking 
AJ = 1 mixings between states of the same I 3 and Y in the 8 and 10, such as A + p. Quite 
simply, these descriptions represent two complete bases of a vector space corresponding to 
a particular class of observable: One basis is organized in such a way as to give one basis 
vector for each observable for a given state, and the other basis is organized according to 
quantum numbers of the spin-flavor symmetry. In such an analysis, an arbitrary amount of 
symmetry breaking can be accommodated. 
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This approach was used to classify all static observables of the literal SU(6) 56 (i.e., 
using only N c = 3) in Ref. with a deeper study of quadrupole moments in Ref. j^l- It has 
been used in the 1/N C expansion several times: for the masses of the 56 ^(|, for charge radii 



and quadrupole moments 



Jul 



baryon multiplet 70 



13 L for the masses and couplings of the orbitally-excited 



and even for the A — > couplings closely related 



to the AN transition magnetic moments [19]. 

This paper is organized as follows: In Sec. ITU we explain why the operator expansion 1/N C 
truncates at a finite order, and how the complete set of operators may be enumerated. We 
compute and tabulate all the matrix elements of all these operators in Sec. IIHI In Sec. IIVI 
we compare our approach to previous ones in the literature [with and without perturbative 
flavor SU(3) breaking], derive new relations, fit to all existing data, and use the results of 
this fit to predict all unmeasured moments. The casual reader uninterested in calculational 
details is encouraged to skip directly to Sec. IIVI We summarize and conclude in Sec. El 



II. OPERATOR BASIS 



Each baryon state belongs to a representation composed of N c color fundamental rep- 
resentations combined into a color singlet. While it is suggestive to think of each such 
fundamental representation being associated with a single current quark, such an identifi- 
cation is not necessary; in general, each fundamental representation merely represents an 
interpolating field whose quantum numbers match those of a single quark in color, spin, and 
flavor — each of these in the fundamental representation of the corresponding group — and 
which together exhaust the whole baryon wave function In general, such a field con- 
sists of not only a current quark, but gluons and sea quark-antiquark pair Fock components 
as well, and indeed may be thought of as a rigorously defined constituent quark. We continue 
to denote such an interpolating field by the simple label "quark." 

An arbitrary baryon bilinear, as appearing in the Hamiltonian for masses, magnetic mo- 
ments, etc., thus carries the quantum numbers of iV c quarks and N c antiquarks. Since 
fundamental (g) antifundamental = adjoint © singlet for all SU groups (in this case, each 
of spin, flavor, and spin-flavor) , each operator that can connect the baryon states may be 
decomposed into pieces transforming as products of 0, 1, 2, . . . , up to N c adjoint represen- 
tations. In terms of spin-flavor SU(2iVV), where Np is the number of light flavors, the 
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operators comprising the adjoint are denned: 



G ia = E^(y®y)^ (2.1) 

where the index a sums over the N c quarks, a 1 are the Pauli spin matrices, and A a are the 
Gell-Mann flavor matrices. Thus, each distinct operator may be written as a monomial in J, 
T, and G of total order n, with 0< n < N c . Such an operator is termed an n-body operator. 

A large subset of operators constructed in this way are redundant or give vanishing 
matrix elements due to group-theoretical constraints. For example, commutators such as 
[J 1 , J 3 '] =ie^ k J k behave exactly as they do for the underlying a and A matrices. Furthermore, 
some combinations of J, T, and G act only on non-symmetric combinations of quarks and 
hence annihilate the ground-state wave functions, while yet other combinations are spin- 
flavor Casimirs and hence give the same value for every state of the representation, making 
them indistinguishable from the identity operator. The operator reduction rule for removing 
all such extra operators was derived for the 56 in Ref. 12(1 . and extended to the 70 in 

n 

Refs. For the present case with Np<3, the rule states: All operator products in which 
two flavor indices are contracted using d~ ab , d abc , or f abc or two spin indices on G's are 
contracted using <5* J or e l ^ k , can be eliminated. 

None of the preceding reasoning depends specifically upon the 1/N C expansion. Such 
1/N C factors arise from two sources: First, an n-body operator appears in an irreducible 
diagram in which n quarks are connected by gluons, requiring a minimum of n— 1 gluons 
exchanged; the 't Hooft scaling a s oc 1/N C then implies an explicit suppression factor 1/iV™ -1 . 
Second, the combinatorics of quarks inside the baryon permits the matrix element of J, T, 
or G to be as large as O(N^) whenever the contributions from the N c quarks add coherently. 
However, if the baryons chosen nevertheless have spins, isospins, and strangeness of 0(N®)— 
as we choose for the spin-^ 8 and spin-| 10 — then the matrix elements of J 1 ' 2,3 , T 1 ' 2,3 , and 
jV s = |(l-2V3T 8 ) are also 0(N°). 

The replacement of T 8 by N s in constructing the operator basis presents a trivial example 
of what, in Ref. fig , is called "operator demotion." Whereas operator reduction rules 
identify linear combinations of operators that give precisely zero to all orders in 1/N C when 



acting upon all states in a baryon spin-flavor representation, operator demotion identifies 
linear combinations of operators whose matrix elements are a higher order in powers of 1/N C 
than those of the component operators, at least for the observed baryon states. Such a 
result can occur since the J, T, and G matrix elements in general contain both leading and 
subleading contributions in N c . 

In summary, a complete accounting of the 1/N C expansion thus requires one to take into 
account the ingredients: i) a complete set of operators under spin-flavor; ii) operator reduc- 
tion rules to remove linearly dependent operators; Hi) a counting of N c factors arising both 
from explicit powers of a s on one hand and coherent contributions due to quark combina- 
torics on the other; and iv) operator demotions to identify operators whose matrix elements 
are linearly dependent at 0(1/N™) for some n but are independent at 0(l/N™ +1 ). 

A full analysis of all states in the baryon multiplets for N c large and finite, as discussed 
above, requires the inclusion of up to iV c -body operators. A parallel analysis carried out 
for N c = 3 states, by the same reasoning, requires only up to 3-body operators. Once the 
physical N c — 3 baryons are identified as states embedded within the N c > 3 multiplets, one 
sees that 4-, 5-, . . ., iV c -body operators do indeed act upon the physical baryons, but give 
results linearly dependent on those of lower-order operators, and therefore may be discarded. 

In the case of magnetic moments for the 56, we have seen that there are 27 independent 
parameters with AJ = 1, AJ 3 = 0, AY = 0, A/ 3 = 0, T odd for N c = 3. The conditions 
AJ = 1 and AJ 3 = require that each operator has a single unsummed spin index i, 
which for definiteness we take to be i = 3. T odd, of course, is the behavior of an angular 
momentum under time reversal; as it turns out, this is accomplished automatically because 
all operators containing the structure constants e 4jfc , f abc , or d abc can be eliminated. The 
conditions AY = 0, A/ 3 = require each unsummed flavor index a to equal 3 or 8. The 
complete set of 27 such operators, including the demotion T 8 — > N s appears in Table |l] For 
those cases in which different order ings of component operators would give different values 
for matrix elements (such as J 2 and G 33 ), the operators are written in a symmetric form. 

In fact, a direct calculation shows that no other operator demotions occur. Consider the 
first 6 operators in Table HI the complete set up to and including 0(N®). If the leading 
coefficient of each of these operators — at O(N^) for G 33 and at 0(N®) for the other 5 — for 
each of the 27 observables is collected into a 6x27 matrix, then one finds that this matrix has 
rank 6: No linear combination of the operators has matrix elements that are only 0(N~ l ). 
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Similarly, no combination of the 17 operators up to O(N x ) is demoted to 0(N C 2 ). 



III. COMPUTING MATRIX ELEMENTS 



We compute the matrix elements of the 27 operators listed in Table U] using only the 
Wigner-Eckart theorem (or its variants) and the total spin-flavor symmetry of the 56 baryon 
states. While the task of computing matrix elements of ra-body operators for states con- 
taining an arbitrarily large number (N c ) of constituents may naively seem to require a large 
amount of group-theoretical technology [e.g., SU(6) 9j symbols], it turns out that all of the 
necessary matrix elements can be reduced to simple SU(2) spin and isospin Clebsch-Gordan 
(CG) coefficients, and nothing worse than an SU(2) 6? symbol needs to be computed. All 
of the necessary tools have been developed in Refs. [12J, [l3j , but we present them here for 
completeness. 

We begin by constructing baryon states in the 56. Since the wave function is completely 
symmetric under exchange of spin and flavor quantum numbers of any two quarks, it follows 
that the collection of all N q quarks of any fixed flavor q must be completely symmetric under 
spin exchange. The spin J q carried by them must therefore have its "stretched" value, N q /2. 

Next, the u quarks and d quarks combine to give a state with I 3 = ^(N u — N d ) = J u — J d . 
The total isospin I is determined by noting that the exchange symmetry property of the 
state under u-d flavor exchange must precisely match that of these quarks' spins, in order 
for the total wave function to be completely symmetric under spin-flavor. It follows that 
J ud = I, where J ud = 3 U + J d . 

In the final step, one simply combines the state of ud total spin J ud = I and isospin 
quantum numbers /, I 3 with the symmetrized strange quarks carrying total spin J s to obtain 
the complete state with spin eigenvalues J, J 3 , where J = J ud + 3 S : 



J J 3 ; II 3 (J u J d J s 



I Js 

73 73 



J 3 JS.JS 1 J u J d 



J 3 

•J 7 



ud 



J r 



Jd jf 



° s ° s I 1 



(3.1) 



where the parentheses denote CG coefficients. Now, in order to compute the matrix elements 
of any particular operator, one need only sandwich it between a bra and ket of the form of 
Eq. (|3.1|) and use the Wigner-Eckart theorem. 
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The basic operators T 3 , N s , J 3 , and J 2 , acting diagonally on baryon states, are easy to 
handle even if they are parts of more complicated operators. On the other hand, 

G 18 = ^(J'-SJl), 

= 1(4-4), (3.2) 

are in general not diagonal and must be handled more carefully. According to Table HI they 
appear in the forms G 38 and G 33 , and as 

j*G t8 = -^WJ 2 -3J-J S ), 

fG l3 = ij.(J u -J d ), (3.3) 
which may be simplified by noting that 

J" J, = -^[(J-J s ) 2 -J 2 -J^ = i(J 2 + ^-I 2 ), 
J-(J„- 3 d ) = (J u + J d + J«)-(J„- 3 d ) = J 2 - 3\ + J a .(J„- J„). (3.4) 

It becomes apparent that only a few nontrivial matrix elements need be computed. Denoting 
the matrix element (J a - J^) as (a(3)(°\ where a and (3 are any two quark flavors, the only 
nontrivial required matrix elements are (J 3 ), (Jd), (J%), (us)(°\ and (ds)(°\ 

Even more simplification is possible, because Eq. (|3.1j) depends on the exchange of u 
and d quarks only through the second CG coefficient, and the factor obtained through this 
exchange is just (— i) J «+ J d-^_ Qf CO urse, the eigenvalues J a , which simply count one-half 
the number of quarks of flavor a in these baryons, remain unchanged from initial to final 
state. The same is true for I 3 = J u — Jd, but the total isospin may change to a value One 
thus finds for an operator O that 

(I'l 3 \0(u <-> d)\II 3 ) = — I 3 \0\ I -I 3 }. (3.5) 

Thus, the only matrix elements that need be computed are (J 3 ), (4), an d (us)(°\ These 
were computed in Ref. |13J] and are reproduced here: 

(4) = ^a.a^^^^,^-!)^^^^^^^^^-^ 



X\/J u (Ju + 1)(2J„ + 1)(2/' + 1)(2J + 1)(2J' + 1)(2J + 1) 

/ 




x . „ J, (3-6) 



(Js) = Sj'^i'iSjljJj'jM^- 1 ) 1 ^ +Js+I yMJs + 1)(2J. + 1)(2J' + 1)(2J + 1) 

/ 




1 J' J 

J 3 -J 3 



(3.7) 



(W (0) = 8j IJ 8 J ,,j 3 5j LJ Jv J 8 Jla jS-l) 1+J+Js ~ J ^ Jd 



x^J u (Ju + 1)(2J„ + 1)J S (J S + 1)(2J S + 1)(2/' + 1)(2/ + 1) 

Jd Ju I \ \ J Js I 



1 V J„ II' J. 



(3.8) 



Note that, in the interest of exhibiting maximal symmetry, the remaining CG coefficients 
have been written as 3j symbols. Despite the fact that a number of their entries are O(N^), 
all the 3j and 6j symbols of interest may be computed using analytic forms appearing in 
the standard text by Edmonds [2]]]. 

The matrix elements for all relevant states are presented in Tables ITTVlIXI Tables ITTklHI 
are lifted directly from Ref. (except for the repair of typos in the S°A matrix elements 
in Table HH. 



IV. RESULTS 

If data existed for all of the 27 observables associated with the magnetic moment sector, 
one would proceed to form a Hamiltonian 

H=-fj,B, (4.1) 

where the operators Oi of Table U enter with unknown dimensionless coefficients q via 

27 

Mz = VoY, C i°" ( 4 - 2 ) 

i=l 

where [Iq is the sole scale in the problem, a mean value of magnetic moments in the multiplet, 
which one expects to be some 0(1) multiple of the nuclear magneton fi^. The 1/N C expan- 
sion provides a reliable effective Hamiltonian if the coefficients q are not larger than 0(1). 
In fact, a number of them may be smaller than 0(1) because certain operators may only 
contribute once SU(3) flavor symmetry is broken. They may also be smaller if dynamical ef- 
fects are present that suppress them below the level predicted by naive l/N c counting. With 
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all 27 observables in hand, one would simply invert the 27x27 matrix whose elements are 
given in Tables IIVHIXI and solve for all q to test this hypothesis. Essentially this procedure 
was carried out for the masses of the 56 in Ref. jiol ]. 



However, the Review of Particle Physics 
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gives unambiguous values for only 10 of the 



observables: magnetic moments of 7 of the 8 octet baryons (fi^o is unknown), the Q~ , and 
the S°A and A + p transition moments. The last of these is extracted from the A — > A/7 

helicity amplitudes Ai and A3 via the standard formula for the Ml amplitude: 

2 2 

Ai + V3As 

^ A+p = 1 * , (4.3) 

V KOlk 

where k ^260 MeV is the photon momentum, from which one finds ua+p = 3.51 ± 0.09 lin- 

r 1 

In addition, we use a recent extraction [23[ /^a++ = 6. 14 ±0.51 [In obtained from an analysis 
of data that has some model dependence, but that respects both gauge invariance and the 
finite A ++ width. We therefore include 11 observables in our analysis. There is also a recent 
experimental determination 24j of fj,&+ = 2-7-i,3 (stat) ± 1.5 (syst) ± 3 (theor) fi N , but due 



to the large theoretical uncertainty we do not use this value in our analysis. 

With only 11 pieces of information to study a system of 27 observables, one must resort to 
using the known quantities to fit the coefficients at the lowest orders of the 1/N C expansion, 
and to use the coefficients so obtained to predict the remaining observables. One may 
then proceed either by i) separating the observables into isoscalar and isovector, as well as 
1 = 2 and 3 isotensor, combinations, or ii) one may employ the electromagnetic nature of 
magnetic moments to construct only operators with a flavor dependence in proportion to 
the quark charges (the "single-photon ansatz"). Since both methods have been employed in 
the literature, we discuss them each in turn. 



A. Analysis in the Isoscalar, Isovector, Isotensor Basis 

The analysis of Ref. j^] (JM) separates operators, and the corresponding combinations of 
magnetic moments, into 1=0 and 1=1 forms. Since the maximal isospin appearing in the 
56 is I (for the A), isotensor combinations with 1=2 and 1=3 are also present: 

1 = 2: - 2// E o +//£-), 0"a++- A*A+ - /^A° + A*A-), (/is*+ - 2/i E *o + /i E «-), 

(/^A+p — ^A°n), (/^E*+£+ — 2/X s *o 2 o + /i£*-£-), 

1 = 3 : (a*a++ - 3yU A + + 3/i A o - Ma-)- (4.4) 
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The JM analysis introduces a leading-order operator X ia , which is equivalent to the 
0(N®) part of G m /N c , and a strange quark spin operator 

f s = hj i -2V3G i8 ). (4.5) 

The JM operator basis then consists of the 6 operators N c X t3 and N s X t3 (I = 1), and J\ 
J l s , N s J l /N c , and N s J l s /N c (I = 0). Since no combinations of these operators have / = 2 or 
3, the combinations in Eqs. (j4.4j) exactly vanish, giving relations 11-16 (JM Table 2). 

In comparison with our Table HJ the choice of JM operators reflects the inclusion of all 
(2) with / = 1 at 0{N]) and O(iV c ), and all (4) with 1= at 0(N°) and 0(N~ 1 ). Since 
there are (as one may readily count) 10 /= and 11 /= 1 magnetic moment combinations 
in the 56, it follows that JM predict 6 isoscalar relations (JM Table 3 S1-S6) that receive 
only 0(N~ 2 ) corrections, and 9 isovector relations (JM Table 3 V1-V7, V8i, and V9i) that 
receive only (^(iV" 1 ) corrections. As expected, we confirm these predictions in our basis. 

The JM analysis makes no use of the electromagnetic behavior of magnetic moments, 
nor of perturbative SU(3) flavor breaking; its analysis can be said to hold in the presence 
of arbitrarily large SU(3) breaking. Thus, operators are organized solely by the 1/N C power 
suppression of their matrix elements. Since 17 operators occur up to and including 0{N^ 1 ) 
while only 11 moment parameters have been measured, it is not yet possible to improve 
upon the numerical analysis of JM using their same scheme. One must therefore impose a 
physically natural flavor structure on the magnetic moment operators, a topic to which we 
next turn. 



B. Analysis Using the Single-Photon Ansatz 

Like all electromagnetic multipole moments, magnetic moments are defined through a 
particular coupling of a bilinear to an on-shell photon. The lowest-order flavor structure of 
the coupling to the photon should therefore be such that each quark couples proportionally 
to its electric charge. In particular, in the limit in which all other sources of SU(3) breaking 
are suppressed, only operators with one unsummed flavor index a may appear, and then 
only in the linear combination (a = 3) + (a = 8)/-\/3 = (a = Q). Specifically, these are the 
forms 

Q = T Q = T 3 + — ^=T 8 , 
V3 
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G lQ = G iZ + -p G iS . (4.6) 
v3 

Implicit in this definition of the quark charge matrix Q is that the quarks assume their 
usual N c = 3 values q u = + | , ^ = q s — — | . In terms of SU (3) flavor hypercharges, 1^ = = | , 
F s = -§. An alternate choice, q u = (N c + 1)/(2JV C ), q d = q s = (~N C + 1)/(2AQ (Y u =Y d = 
1/N C , Y s — — 1+ 1/N C ), has the convenient property that all hadrons then have the same 
electric charges and hypercharges for arbitrary iV c as they do for N c = 3. Moreover, with this 
choice the chiral anomalies of the Standard Model (with N c colors) automatically cancel. 
However, one is also faced with the mysterious prospect of electromagnetic charges dependent 
upon the number N c of QCD charges. More significantly, the quantization condition of 
the Wess-Zumino term permits only baryon SU(3) representations containing states with 
hypercharge Y = N c /2> j25|; if such states have 0(N®) strange quarks, then the ^-dependent 
choice (Y u = Y l i = l/N c ) is disallowed. For the remainder of this paper, we assume the usual 
A^-independent quark charges. 

The only operators occurring in the single-photon ansatz with no other SU(3) breaking 
(cf. Table UJ) are 

0\ = G 3Q , 2 = ^QJ\ 6 3 = ~{J 2 ,G 3Q }, 4 = ^-J l G tQ J 3 . (4.7) 

For any value of N c it turns out that the combination (9 3 = ((%— O4) vanishes for all diagonal 
moments and survives only for transitions. Since the only transition moment measured at 
present is fj,&+ p , using O3 rather than O3 provides a more incisive test of the expansion when 
fitting to current data. 

In addition, one may perturbatively break SU(3) symmetry by including effects due to 
a finite strange quark mass or spin. We incorporate such effects by including a parameter 
e that indicates each instance of breaking of the SU(3) symmetry. At first blush, one may 
suppose that it is proportional to m s , e ~ 0.3 ~ 1/N C , but as we see below such a rigid 
identification is not necessary. The list of additional operators with matrix elements up to 
0(e x JVj) reads 

eO s =eq,Jl eOs = jj-N a GW, e(D 7 =-^QJ 3 s . (4.8) 

These forms are obtained by the simple expedient of inserting sources of SU(3) breaking 
along the strangeness direction into the operators of Eq. (|4.7jl . and retaining only those with 
matrix elements up to O^N®). The possible substitutions are J 1 — > eJ l s , Q/N c — * eq s N s /N c , 

11 



qiQ eq s J l s , or multiplication of an operator by eN s /N c ; each of the last three of these 
replacements also costs a power of 1/N C . In terms of this basis, the analysis of Ref. 
(LMRW) uses the operators O x [LMRW Eq. (17)], 2 [LMRW Eq. (21)], the combination 



(4.9) 



N c 

from chiral loop diagrams [LMRW Eq. (25)], and the operators m s 5 , m s 6 , m s 7 from 
counterterms to the loop calculation [LMRW Eq. (29)], 5 independent operators in all. Note 
the characteristic nonanalytic m s behavior in Eq. ([4.9)1 . which suggests that the appropriate 
SU(3) expansion parameter e might properly scale as ^iV" 1 / 2 rather than ~ A" 1 ; we address 
this issue below. 

At the next order, 0(e 1 N~ 1 ), one finds the 6 operators 

e0 8 = eq s ^J 3 , £0 9 = s^QJ 3 , sO 10 = ~{J- J s , G 3Q }, 

eOn = ^-/ s G jQ J 3 , e0 12 = ~{J j G> Q , J 3 }- (4-10) 

Beyond this collection, the next operators have matrix elements of 0(e 2 N~ 1 ) and 0(e 1 N^ 2 ), 
but significantly there are none of 0{e 2 N®) or 0(e°N~ 2 ). As a consequence, whether one 
takes the SU(3)-breaking parameter e ~ N^ 1 ^ 2 or ~ iV" 1 , the series expansion truncates 
consistently after the inclusion of either the set Ci, . . . ,e0 7 (complete to combined orders 
e 1 ^ and e ^ 1 ) or the larger set O u . . . ,eO n [complete to 0(e l N- 1 )}. 

In this notation, the operators used in Ref. Ill (DDJM) consist of [DDJM Eq. (2.9)] 

01, 2 , 3 , £05, £06, £07, £08, (4.11) 

and the operator e 2 q s N s Jl/N c with a coefficient fixed relative to those of Oq and 07, 7 
independent operators in all, but a somewhat different set than Oi, . . . ,eOj. Note that 
DDJM does not assign particular powers of e; DDJM also recognizes the presence of 0{w}J 2 ) 
loop corrections, so that statements regarding the meaning of e still apply. 



C. A Global Fit 

As we have seen, assigning a particular numerical value to the SU(3)-breaking parameter e 
can be problematic, owing to the existence of contributions nonanalytic in m s . Fortunately, 
we have also seen that regardless of whether one takes e ~ A" 1 / 2 or A" 1 , the expansion 
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truncates consistently after the 7 operators in Eqs. (|4.7|) . ()4.8|) or after the 12 operators 
including Eqs. ()4.10|) . Powers of e may simply be left in the coefficients, in which case the 
size of SU(3) breaking for each operator may be judged directly from a fit to data. 

Unfortunately, with only 11 measured magnetic moment parameters, only a fit to the first 
7 operators is possible at present. We therefore perform a least-squares fit to the coefficients 
in the expansion 

7 

i=l 

We choose the scale /xq by reasoning that the best known value among the moments is fi p , 
and that there is only one operator, Ox, at leading order (vV*), whose value for the proton 
is (iV c + 3)/12 = 1/2 for N c = 3. Therefore, a natural choice that makes the sole leading- 
order coefficient d\ of order unity is to set /io = 2/x p [In alternate choices one may average 
over several measured magnetic moments, but this merely renormalizes all d{ by the same 
0(1) multiple]. Since the expansion is truncated by ignoring effects of 0(e 1 N~ 1 ), one must 
include in addition to the statistical uncertainty of each magnetic moment a "theoretical 
uncertainty" of order [iqe/N c . In fact, the x 2 /d.o.f. obtained from a theoretical uncertainty 
of H P /Nl ~ 0.3 /xjv; f° r example, is only about 0.13, suggesting that the naive theoretical 
uncertainty is a gross overestimate. We find empirically that choosing it to have a value 
about n P /Nl ~ 0.1 fiN gives a x 2 /d.o.f. = 1.05, meaning that the fit is as good as one might 
hope for even if all 0{e l N~ 1 ) effects are suppressed, and uncertainties are effectively only 
0(e 2 N~ 1 ). This result far supersedes that expected from a naive 1/N C expansion. 

The fit values for the coefficients are given in Table [X] One immediately notes that no 
coefficients are larger than 0(1); had any of them been so, one would conclude that the 1/N C 
expansion fails. But in fact, d\ and d% are of 0(1), while c?2 and d± are actually consistent 
with zero. The SU(3)-breaking coefficients g?4,5,6, do indeed have central values about 1/3 
or less, but only d^ is statistically different from zero. Such a pattern of suppression beyond 
naive 1/N C counting has in fact been seen before, in the orbitally excited baryons [15I lid Av\ . 
Moreover, a hint of this effect is visible in the results of DDJM Table VI [although their 
operator basis, and especially their treatment of SU(3) breaking, is rather different]. We do 
not understand the source of this suppression beyond that expected from l/N c counting, 
and find it to be the most intriguing feature of our analysis. 

One may also use the fit values for di to predict all of the remaining 16 magnetic moments 
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to within the theoretical uncertainty; the results are presented in Table IXII Note that the 
recent //a+ measurement easily agrees with our prediction. An important point that may 
not be obvious from this compilation is that, with only 7 operators included, there remain 
20 relations among the magnetic moments. Among these are vanishing of the 1 = 2 and 
3 combinations in Eq. (J4.4j) : The leading-order operators in Eq. (J4.7|) contain only 1 = 
and 1 pieces, while inserting SU(3) breaking along the strangeness direction induces only 
1 = corrections. The combinations in Eqs. ()4.4|) receive contributions only from tiny 
isospin-breaking effects due to either {m u — m d ) or loop diagrams containing an additional 
photon [0(a/47r)]. Furthermore, /ia« vanishes for all 7 operators when N c = 3, and receives 
contributions only from SU(3) breaking not solely in the strangeness direction (since A 
contains no s quarks), which induces an additional m u ^/rn s suppression factor. A similar 
statement holds for all nonstrange observables: None of them receive a contribution from 
any operators beyond 0±, . . . , O4; indeed, the famous SU(6) relation /i n = — 2/i p /3 for N c = 3 
receives a contribution only from the anomalously suppressed operator 02. 

Predictions of the diagonal moments appearing in LMRW Table I agree well with the 
results of Table IXl) but ours have smaller uncertainties due to the larger number of operators 
and the treatment of subleading effects as described above. 

Since the next order of the expansion contains 12 operators while 1 1 observables are well 
known, it is tempting to suppose that just one more moment measured — say, an improvement 
on — w iH permit such a fit. However, the I = 3 and 1=2 relations among the A's 
combined with the result n^o = 0, as satisfied by all operators in our list [and any others 
breaking SU(3) solely in the strangeness direction], predict that the A magnetic moments are 
exactly proportional to electric charge for N c = 3: /xa++ = 2/ia+ = —2[j,&-. Moreover, there 
is precisely one relation satisfied by the first 12 operators among the measured moments: 

1 3 

f*n ~ + £*£-) - ~ v^/Usoa + A*3° = 0(e 2 N~ 1 ), (4.13) 

which has a numerical value of 0.22 ± 0.14 n^, the uncertainty being completely dominated 
by that of p,y,o\. In fact, Eq. ()4.13|) was originally derived in heavy baryon chiral perturbation 
theory |26j], where it was found to have no 0(ml^ 2 ), 0(m s ), or 0(m s \nm s ) corrections — i.e., 
no 0(e l ) corrections in the current formalism. Converting Eq. ()4.13j) into a scale-invariant 
result by dividing by the average of the same expression with all negative values turned to 
positive ones [giving an O(N^) combination], one obtains 0.057 ± 0.036, in good agreement 
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with expected magnitude e 2 N~ 2 . A better measurement of /i£°A ; certainly within current 
experimental means, would decisively test the expansion at this order. 

It follows that a measurement of at least two moments with nonzero strangeness is re- 
quired to perform a fit to the 12 operators at 0(e 1 N~ 1 ), and to determine whether the 
effects at this order are truly suppressed, as the x 2 /d.o.f. suggests. The S*A, S*S, and S*S 
transitions are natural candidates, since the associated radiative decays are presumably be- 
ing recorded (although not yet studied) at Jefferson Lab, as well as other facilities. To date, 
the decay S*° — > A7 has been seen in precisely one event (27 1; the opportunity to improve 
on this meager set clearly exists. 

Finally, we note in passing that once the unmeasured moments are included, a number 
of relations with only 0(e 2 N~ 1 ) corrections, in addition to Eqs. (|4.4j) . (|4.13|) . and = 0, 
remain (7, to be precise). A particularly elegant example is = 



V. CONCLUSIONS 



We have developed a basis of operators representing every possible observable pattern 
of magnetic moments for the ground-state spin- 1/2 and spin-3/2 baryon multiplets, and 
organized them according to the counting of 1/N C factors. We have furthermore computed 
the group-theoretical parts of all of these operators, thus producing a complete effective 
Hamiltonian for magnetic moments. Our analysis of this operator expansion examined the 
consequences both in the case of arbitrarily large SU(3) breaking and perturbative SU(3) 
breaking (in powers of a parameter e) beyond that produced by coupling quarks to photons 
in proportion to their electric charges. In both cases we have compared to previous results 
and showed how this work extends earlier analyses. 

In particular, we have found in the case of nonperturbative SU(3) breaking that the mea- 
surement of several additional magnetic moments is necessary to improve numerically upon 
previous analysis [i.e., from relative order N~ 2 to A"^ 3 ]. However, in the more physically 
meaningful case of perturbative SU(3) breaking, the series may be truncated consistently 
after 7 operators (including up to orders e^AT? and e°N~ 1 ) or after 12 operators [up to 
0(e 1 N~ 1 )]. Since 11 observables are currently well measured, we presented results of a fit 
to 7 operators, and found not only that several of the effective Hamiltonian coefficients are 
smaller than expected, but also that a good fit can be obtained if the terms neglected are 
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actually 1/N C smaller than naively expected. 

A number of relations among the magnetic moments survive the expansion to 12 oper- 
ators. After enumerating a number of them [e.g., Eq. (|4.13|) ]. we suggested that the most 
probative tests of the 1/N C expansion: an improved measurement of the S°A transition and 
observation of £*A, and S*S transitions, should lie with current experimental means. 
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TABLE I: The 27 linearly independent operators contributing to the magnetic moments of the 
spin-i and spin-| ground-state baryons, organized according to the leading N c counting of their 
matrix elements. 
0{Nl) G 33 

0(iV c °) J 3 , G 38 , j^T 3 G 33 , -±- c N s G 33 , ^{J^G i3 ,G 33 } 
O(iV-i) ^-T 3 J 3 , ^-N S J 3 , ±T 3 G 38 , ±N S G 38 , ^1{J^,G 33 }, ^{T 3 ?G 33 , ^N*G 33 , 

0(iV- 2 ) ^J 2 J 3 , ^iV s 2 J 3 , ^(T3) 2 J3, ^T3JV S J3, ^{ji,G 3S }, ^(T 3 ) 2 G 3 ^iV 2 G 38 , 
■^T 3 N S G 38 , ^J*G iS J 3 , ^l{J*G i8 ,G 38 } 
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TABLE II: Matrix elements of the operators N uAs [whence = (aa)(°) = (N a /2)(N a /2 + 1)] 
and the rank-0 tensors (a(3)^ with a / (3. Since spin is unchanged by these operators, the matrix 
elements vanish for all off-diagonal transitions except S°A; in that case, the only nonvanishing 
entries are (us)(°) = -(ds)^ = -§ y/(N c - 1)(A C + 3). 
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TABLE III: Matrix elements of the operators J^, J^, and in the state of maximal J 
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TABLE IV: Matrix elements of the 1-, 2-, and 3-body operators corresponding to magnetic mo- 
ments. 
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TABLE V: First continuation of above table. 
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TABLE VI: Second continuation of above table. 
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S° 








Q 
O 





A 








V3 
8 





S°A 





-^^-1)^ + 3) 








E~ 


l 

2 




8 


-n(^c + i) 


H° 


1 

2 




3^ 
16 






_ 1 

z 


4^ 

4o 


_3yj? 

1 ft 
10 


144 


A+p 





^V(^e-l)(iV C + 5) 





i#c-l)(^c + 5) 


A°n 





-3=^-1)^ + 5) 





-^(^-1)^ + 5) 


S*°A 





^(^-1)^ + 3) 
















3 [3 

A \ 9 












3 [3 

4 V 2 


12^(^ + 1) 


£*-£- 







3 /3 

4 V 2 











3 /3 
IV 2 


1 N 

36V2 c 









3 [Z 

4 V 2 


36^ 
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TABLE VII: Third continuation of above table. 



State 


((T 3 ) 2 G 38 > 


(A/ 2 G 33 > 


<7V 2 G 38 ) 


(T 3 N S G 33 ) 


A++ 


9^ 
16 











A+ 


id 











A 


ID 











A" 


16 











£*+ 





\{N C + 1) 





h(N c + 1) 


S*° 














s *- 





-l(N c + l) 





g(^c + l) 




^3 
16 


lN c 








16 


-$N e 




T2 N c 


n- 








2 





p 


1 

WV3 











n 


1 

16^/3 











£+ 


1 




1 


n(JVc + l) 


S° 








1 

2\/3 





A 








1_ 





S°A 





-^V(^c-l)(iV c + 3) 








E" 


1 

2V3 


-n(^c + i) 


1 




E° 


Vs 

16 




-Vs 






VI 

ID 


y 


-V3 


oD 


A + p 














A°n 














S*°A 





-W(#c-l)(JVc + 3) 

DV ^ 
















1 
V6 







1 


12 U(^ +1 ) 


1 

V6 


1^(^+1) 


£*-£- 


1 

v& 


12 y^c+l) 


1 






1 




2V2 






1 




2V2 


973^ 




4^6 


V3 
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TABLE VIII: Fourth continuation of above table. 



State 


(T 3 A S G 38 ) 




((J^V 3 ) 


^({J^G 33 }) 


A++ 





U N * + 2 ) 


15^ 
16 


•^(A c + 2) 


A+ 





iK^c + 2) 


15^ 
id 


#(iV c + 2) 


A 





-^(N c + 2) 


Id 


-!§(N c + 2) 


A~ 





-^(iV c + 2) 


16 


_3^3 (A r c + 2 ) 


£*+ 





i(iv c +i) 








S*° 














£*- 





-i(iv c +i) 








■^*o 


4 


-^-A„ 

16 c 


15^ 
16 


5 N 




V3 
4 


5. /y 

16 c 


15^ 
16 


5 A 

on, /q C 


n- 








15\/3 
8 





p 





ili(^ c + 2) 


V3 
16 


32^ + 2) 


n 





-n(^c + 2) 


16 


-^(A c + 2) 

32vo 


£+ 


1 

o /q 


n(^c + l) 


8 




S° 








V3 
8 





A 








V3 
8 





S°A 





--Ly(jV c -l)(jV c + 3) 








E- 


1_ 


~U N c + 1 ) 


8 


-1673^ + ^ 

IDV O 


H° 


4 




3^ 
16 


on /q-^C 




4 




16 


1 A 

32\/3 c 


A+p 











11^(^-1)^ + 5) 


A°n 











WiVw=-i)(tf c +5) 


S*°A 











16 ^V(^ C l)(JVc + 3) 


















1 










£*-£- 


1 

VE 








32V6^ + 1 ) 




1 








L_ _/\r 

6^ c 




1 








_J_/V 

ev^ c 
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TABLE IX: Fifth continuation of above table. 



State 


±({J<G i8 ,G 38 }) 


±({J*G i3 ,G 33 }> 


±({J*G i3 ,G 38 }) 


A++ 


15 

32 


lio(iV c + 2) 2 




A+ 


15 

QO 




#(iV c + 2) 


A 


15 

32 


m( N c + 2 ) 2 


-^|(Ar c + 2) 


A- 


15 

32 


lio(^ + 2) 2 


_^3 (A r c + 2) 


£*+ 





Tii(^ + 1 ) 2 





S*° 











T,*- 












15 
32 


-JL i\T 2 
288 c 


5 N 

3^/3 C 


£*- 


15 
32 


-JL jV 2 

288 c 


5 iV 


n- 


15 

8 








V 


1 

32 




^7f(^Ve + 2) 


n 


J_ 


UNc + 2) 2 

yo v ' 


_ 1 (AT c + 2) 

32V3 v y 


£+ 


1 

8 




1673^ + 1 ) 

ID V «5 


S° 


1 

8 


±(N c -l)(N c + 3) 





A 


1 

8 


±(N c -l)(N c + 3) 





S°A 











E" 


1 

8 




-1673^ + ^ 
lov 


H° 


9 

32 


J_ i\r 2 

864 c 






9 

32 


J- AT 2 

864 c 


1 jV 

32\/3 c 


A+p 





-J 7 =(N c + 2W(N c -l)(N c + 5) 





A°n 





-j^(N c + 2)^(N c -l)(N c + 5) 





E*°A 










E*°E° 


1 

8 a/2 








E*+E+ 


1 

8i/2 


384^ + 1 ) 2 




E*-E- 


1 

SV2 


384 7 v ^(^ + 1 ) 2 


321/6^ + i ) 


~*°~° 


\_ 

2V2 


1 AT 2 
108^ c 


1 N 
12VE c 




1 

2^2 


1 N 2 

108 V2 c 


12^ 
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TABLE X: Best fit values for the coefficients in the expansion Eq. I|4.12|) . 



di 


= +0.987 ± 0.038 


d 2 


= +0.076 ± 0.092 


d 3 


= +1.385 ± 0.258 di = +0.059 ± 0.255 


d 5 


= -0.348 ±0.114 


d c> 


= -0.140 + 0.108 


d 7 


= +0.126 + 0.108 



TABLE XI: Best fit values for the 16 unknown magnetic moments in units of /xjy. 
// A+ =+3.04 ±0.13 Mao = +0.00 + 0.10 fj, A - = -3.04 ± 0.13 /i S .+ = +3.35 ± 0.13 
= +0.32 + 0.11 /i E ._ = -2.70 ±0.13 ^~*o = +0.64 ± 0.11 ,u 3 *- = -2.36 ± 0.14 
= +0.77 ±0.10 MAOn = +3.51 + 0.11 me.oa = +2.93 + 0.11 /U S *o s o = +1.39 ± 0.11 
M S , +E + = +2.97 + 0.11 pe-e- = -0.19 ±0.11 // a .o s o = +2.96 + 0.12 // a .- s - = -0.19 ± 0.11 
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